Low Energy Properties of the (n,n) Carbon Nanotubes 
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According to band theory, an ideal undoped (n,n) carbon nanotube is metallic. We show that the 
electron-electron interaction causes it to become Mott insulating with a spin gap. More interestingly, 
upon doping it develops superconducting fluctuations. 
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Carbon nano-structures such as Cqq |lj and nanotubes 
H have attracted considerable interest recently. The lat- 
ter is a graphite sheet wrapped into a cylinder form. A 
a pair of integers (n, to) specifies the wrapping. Starting 
from a graphite sheet with the primitive lattice vectors 
a, b making an angle of 60 degrees, the (n, to) tube is a 
cylinder with axis running perpendicular to na + mb, so 
that atoms separated by na + mb are wrapped onto each 
other. 

Considerable efforts have gone into studying the band- 
structure of carbon nanotubes || . The purpose of this 
work is to address the effects of electron-electron interac- 
tion on the low energy properties of them. In one-electron 
tight-binding description where one retains a single it 
orbital per atom and keeps only the nearest neighbor 
hopping, all n ^ to tubes are band-insulators with gaps 
generally scaling inversely with the radius ||| . For these 
tubes a sufficiently weak electron-electron interaction is 
not expected to change the low energy properties qual- 
itatively. The same can not be said of the (n,n) tubes, 
which are band metals. To be precise, for the latter two 
out of An bands intersect the Fermi level to form two 
Dirac points (Fig.l). Due to the low dimensionality and 
the presence of gapless excitations, the effects of electron- 
electron interaction must be examined more carefully. 



+ 





FIG. 1. The low-energy band structures of (n,n) tubes. 
"L" and "R" labels the Dirac points, "+", "-" labels the right 
and left movers. The dashed line denotes the Fermi level in 
the doped case. 

In this paper we perform perturbative renormaliza- 



tion group (RG) calculations to analyze the asymptotic 
low energy behavior of the (n,n) tubes (40. The mi- 
croscopic Hamiltonian is consisted of a nearest-neighbor 
tight-binding model for the 7r orbital on each carbon atom 
(the bonding topology is illustrated in Fig. 2) and a 
Hubbard U plus a nearest neighbor V for electron cor- 
relations. To build up an effective model for low tem- 
perature, we first discard all bands that do not intersect 
the Fermi level. Second, we regard the two remaining 
bands as linearly dispersing, i.e., we ignore the energy 
dependence of the band velocity. Due to these approx- 
imations, a upper energy cutoff E c has to be imposed 
on our subsequent discussions ||. In the Hilbert space 
of the two bands, the original interactions U and V give 
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rise to twelve independent scattering amplitudes g^ cd 
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FIG. 2. The bonding structure of (n,n) tubes. 

The effective Hamiltonian is given by H = Hk + Hi, 
where 
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dx ^hA x )^]b*> WVW' (x)yj Ma (x). (l) 



In the second line of Eq.(|l|) the sum is taken over 
i,j,k,l — R,L, a,b,c,d — ±, and a, a' =],[■ The op- 
erator ip^ i _ rT (k) = J dxe lkx ! ip i±a (x) creates an right/left 
moving electron with momentum ki + k and with spin 
a. Here ki is the momenta associated with the right (R) 
or left (L) Dirac points (Fig.l). The only nonzero g's in 
Eq. (m) are those whose lower and upper indices have the 
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form: (RRRR), (RLRL), (RLLR), and (+ + ++), (+ + 

), (H 1 — ), (H h), plus those generated by + «-> — 

and/or R «-> L. We abbreviate the twelve independent 
coupling constant as g^, where i = 1,2,4; j = 1,2,3,4 

§■ 

We then perform a one-loop RG calculation, which is 
a straightforward generalization of the one-band calcula- 
tion The only complication is that instead of four, 
there are twelve independent scattering amplitudes. This 
calculation yields: 
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functions are chosen because of their logarithmic diver- 
gence in the absence of interaction. To the lowest order 
in g's, the results are 
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Here Xo = 5^7 
cept that Nqdwi 
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ln( — ), and all N a are unity ex- 
2. Finally the expo- 
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I. The undoped case. The bare values of the cou- 
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pling constants are g\ — gl = g\ = 

g\ =g\ = g\ = g\ = u, gl = g\ = gl = g% = u 

solve Eq.(|J) numerically for V/U < 1 and find that in all 
cases the absolute values of all coupling constants even- 
tually diverge. As to the susceptibilities, we find that 
for jj < {jf)c ~ 0.8 the most divergent susceptibility 
is Xsdwi, and for (jj) c < jj < 1 the most divergent 
susceptibility is XCDW2- 

Ia)The SDW1 Phase: For ^ < (^) c we perform 
a mean-field theory using the renormalized Hamiltonian. 
Guided by the susceptibility result, we introduce M R = < 
O r (x) > and M L =< O l {x) >, where 

Or(x) = ^[ipR+^tpR-tix) - ^ +l (x)i/jR- L (x)) 
1 



3434 + 23434 - 23232 + 3132 + 3231 +3x3!, (2) L (x) = -[^(x)tJj l ^(x) - t/j+^x)^^)} (5) 



Plus = = = 0- In the above x ee ^H§), 
where E c and E are initial and running energy cutoffs 
respectively. 

To draw the implications from the RG equations we 
compute correlation functions Xq(^) = / dxdte luJt < 
T[O a (x,t)0+ (0,0)] > for the following O a : 

Ocdwi(x) = ^2i/j R+a (x)ifjR-*(x), 

(J 

S dwi{x) = ]^[ip R+] {x)ij R ^{x) - ip R+L (x)ip R ^i(x)] 

Ossi{x) = -^[^fl +T (aj)^_|(a:) ~ ^R+ii^R-^x)}, 
Otsi(x) = ip R+ ^(x)i/; R _^(x), 

C DW2(x) = ^i) R+a (x)i> L - a (x), 
a 

S dw2(x) = ]^[ip R+] {x)i; L ^{x) - \j) R+i (x)ijj L -i{x)), 
Oss2(x) = -j=[ip R+r (x)'ip f L _ l (x) - ipR +L (x)^l_ t (x)], 
T S2(x) = W R+ i(x)tp' L _^(x). (3) 

plus those generated by R «-> L and/or + «-> — . In the 
above CDW, SDW, SS and TS stand for charge density 
wave, spin density wave, singlet superconductivity and 
triplet superconductivity respectively. These correlation 



as the order parameters. The mean-field Hamiltonian is 
the factorized form 

H = H K - g\ J dx[0 R R + 0+0+ + L L + 0+0+] 

- 2g\ J dx[0+0 R + 0+O l ] - 2gl[0 L R + O+O+] 

-2g\ J dx{0+0 R + 0+0 L ]. (6) 

In the range of V/U considered here the renormalized val- 
ues of all the g's appearing in Eq.(||) are positive. Con- 
sequently, the mean-field solution predicts M R = M R = 
AIi = A/£ = Mq. It is important to note that the term 
proportional to gl demands both M R and to be real 
and the term proportional to g\ requires them to have 
the same sign. Consequently, the mean-field theory com- 
pletely fixes the U(l) phases of the order parameters M R 
and Ml- The only global degree of freedom left is the 
SU(2) rotation of the order parameter away from the spin 
z-direction. The Goldstone mode associated with the lat- 
ter governs the low energy physics of the SDW1 phase. 
Indeed, let us define 

Qaa'{x) =< lp R+a (x)lp R - a i(x) - -^-W R+T {x)tl)R- T [x) > 



^ L+T {x)^L-r{x) > +h.C. (7) 
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The repeated indices imply summation. The mean field 
solution corresponds to Q aa > = Q%„i = 2aAloS aa >. A 
smooth twist in the direction of the order parameter cor- 
responds to Qoo' — [U + Q°U} aa > , where U is a smooth 
space-time dependent SU(2) matrix. After a proper 
rescaling of the space and time coordinates, the effective 
action governing the dynamics of such twists is 



S a = | J dtdxTr[{d^Q) 2 }. 



Unlike the half-filled one-band Hubbard model, the ef- 
fective action does not contain a topological term. This 
is because there are two Dirac points instead of one. In- 
deed, it can be shown that each Dirac point contributes a 
topological term with coefficient 7r Q . Thus when added 
together, their effects cancel! Q The 1+1 dimensional 
non- linear a model given by Eq.(||) is always disordered 
f[of , which indicates a gap in the spin spectrum jll| . 

Ib) The CDW2 Phase: For 1 > £ > (^) c , the 
most divergent susceptibility is xcdw 2 - Therefore we 
introduce Ai =< 0\{x) > and A 2 =< 2 (x) > where 
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2 (x) =^2i}l +a (x)ip R - a (x), 



(9) 



as order parameters. The mean-field Hamiltonian is the 
factorized form of 



H= H 



K 



9l) 



--(2 ffl 1 - 52 2 )[0+0 1 



dx\O x 2 

o+o 2 \. 



ofot) 



(10) 



For the range of V/U considered here the renormalized 
(2(7^ — g!) and (2(7^ — c/^) are positive and negative respec- 
tively. In the mean-field theory the term proportional to 
(2gf — g 2 ) sets the sum of the phases of Ai and A 2 to tt, 



i.e., Ai 



-A 2 . The mean-field solution is not unique, 



there being a family of equivalent mean- field solutions re- 
lated by the Ai — * e l8 Ai and A 2 — > e~ l8 Ai transforma- 
tion. This continuous degeneracy reflects the symmetry 
of Eq. (^0|) under this transformation. Under such circum- 
stance there will be a Goldstone mode. To determine the 
effective action for the latter, we let Ai(a;, t) = Aoe 10 ^'^ 
and A 2 (x,t) = —Aqc^ 18 ^'^ in the mean-field Hamilto- 
nian and integrate out the electronic degrees of freedom 
||. (In the above 9(x,t) is a smooth space-time func- 
tion). To determine whether or not this mode is charged, 
we impose a background electric field E. The details of 
such calculation will be reported elsewhere ||, but the 
following is the answer. After proper rescaling of x and 
t the effective action for 6 is 



S -* 
2 



J dtdxty*^ 



(11) 



where <f> = e ,s '. The lack of electric field dependence in 
Eq.(O) indicates that the Goldstone mode is neutral. 



II. The doped case. In the doped case g\ = g\ — 
#2 = 0. Integrating Eq.(||) numerically for V/U < 1, 
we see that the absolute values of all coupling con- 
stants again eventually diverge. In this case there are 
three phases. For y < {jj) c i ~ 0.55 the most di- 
vergent susceptibility is the superconducting XSS2^- For 
{jj)ci < jj < (jj)c 2 ~ 0.65 the most divergent suscep- 
tibility is xsdwi- Finally for {jj) c2 < jj < 1 the most 



(8) divergent susceptibility is xcdwi 



Ha) The SS2 phase For ^ < (^)d we introduce 
Ai =< bi(x) > and A 2 =< b 2 (x) > as the order param- 
eters where 

b+(x) = -j=l^(x)ipl u (x) - 4_ ; (z)V>l +T (z)]. (12) 
The mean-field Hamiltonian is the factorized form of [fL2l , 



H= H K + {g\ +g\) J dx[bt{x)b x {x) + b+(x)b 2 (x) 
+{9l + 9l) I dx[b+(x)b 2 (x)+b+(x)b 1 (x)}. (13) 



For the range of V/U considered here, after sufficient 
steps of renormalization g\ 4- g\ becomes negative. Mean- 
while, g\ +g\ stays positive. In the mean-field theory the 
term proportional to (gf + g\) sets Ai = — A 2 . There is 
also a continuous family of mean-field solutions generated 
by the Ai j2 — ► e^Ai j2 transformation. This degeneracy 
reflects the global U(l) invariance of Eq.(fl3"|). The Gold- 
stone mode associated with <fi dominates the low energy 
physics of the SS2 phase. If we let A{x,t) = A e l9 ^^\ 
after proper rescaling of x and t the effective action for 
the Goldstone mode is 



S' 



| j dtdxW^f 



2A„ 



(14) 



Here 4> = e and A^ is the external gauge field. The 
factor of two in front of reflects the fact that the 
Cooper pair is doubly charged. 

Now we address the effect of impurities on supercon- 
ductivity. Nanotubes often have impurities. For a single 
impurity the following terms are added to the Hamilto- 



H imp = E <fc(0)^6«r(0) + M- (15) 

In the above "0" is the position of the impurity. The 
one-loop RG equations for m°' are: 



LR 



d -^ = \{gl+gl-2g\-2gl)u R + R , 



(16) 



plus six other equations obtained by R «-> L and + «-> — . 
All other u's do not renormalize. For positive bare u~[ R 
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and u RR , u~[ R eventually grows upon renormalization, 
while u RR eventually shrinks. Nominally, we would drop 
the latter and keep only the former. However, within 
the range of coupling constants where we can trust our 
perturbative analyses, u2 R { u r%) is onr y amplified (sup- 
pressed) by roughly a factor of two. For this reason, in 
the following we shall analyze the effects of all impurity 
scattering amplitudes on Cooper pairs. 

We let an initial Cooper pair state £> + (0)|0 >= 
(b+(Q) - b+(0))\0 > be scattered by vaious rnipu- 
scattering channels. The following are the re- 
sults. 1) The terms introduced by u~^ R , u£r, u~^~[, u rl 
annihilate the pair. 2) The terms introduced by 

U RRi U RR.> U LL> U LL> U lli U LRi U RL2^RL break _ the P air - 

3) The terms introduced by u RR , , u LL scatter 
it. Therefore a large amount of impurities can destroy 
superconductivity through 2), and produce Copper pair 
localization through 3). 

lib) The SDWl phase For (^) cl < £ < (^) c2 
our order parameters are the same as those defined in 
Eq.(||). The mean-field Hamiltonian is the factorized 
form of Eq.(||) except that gf and gl are set to zero. 
For the range of V/U considered here the renormalized 
values of g\ and g\ are both positive. Consequently, the 
term proportional to g\ requires Mr = Ml- In this case, 
in addition to the global SU(2) freedom associated with 
the spin rotation, there remains a global U(l) freedom, 
i.e., Mr.l — > e^MR.L. The effective action for the spin 
Goldstone mode is the same as that in Eq.(||), hence the 
spin excitations remain gapped. The effective action for 
the charge Goldstone mode is that of Eq.([TT|) excepts 
that an additional term arising from the chiral anomaly 
should be added. 



Sn 



S c hi 



dxdtE8(x,t) (17) 



The last term reflects the fact that the U(l) Goldstone 
mode is a charged mode. The physical meaning of this 
mode is the sliding degree of freedom of the spin density 
wave. 

lie) The CDW1 phase For (£) c2 < j. < i we 
introduce =< Or{x) > and =< Ol{x) > as 
order parameters. In the above 

R {x) =^2ip R+a .(x)i/jR^ cr (x) 

(7 

L (x)='£il>l +a (x)iP L - lT (x). (18) 

a 

The mean-field Hamiltonian is the factorized form of 
H = H K + \{2g\ - g\) J dx[0 R 0+ + O l O+). (19) 

In the range of V/U considered here 2g\—g\ < 0. Conse- 
quently, mean-field theory gives Ar = Al with a global 
U(l) freedom, namely Ar.l — > e l9 Ar^. The Goldstone 
mode associated with this U(l) freedom is again governed 



by an action of the form Eq. ( |17| ) . In this case <fr describes 
the sliding mode of the charge density wave. 

In this paper we use an on-site and nearest-neighbor 
interaction to model the screened Coulomb interaction. 
The effects of an unscreened Coulomb potential remain 
to be studied. In addition, we have not yet stud- 
ied the nature of various phase transitions as a func- 
tion of V/U. Finally, our recursion relation (Eq.(g)) 
has been obtained previously in studying various ver- 
sion of two-band models . The relation between our 
model and the latter can be established via the identifi- 
cation that (R,+) <-> (A,R), (R, -) <-► (B,L),(L,+) <-> 
(B,R),(L,—) *-* (A,L), where A and B labels the two 
bands. 
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